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1. Introduction 


Matrix models a proven to be a powerful tool in the study of various mathematical 
and physical problems, such as random geometry and enumeration of graphs, strings and 
two dimensional quantum gravity, chaos and mesoscopic systems, statistical mechanics of 
spins on random lattices etc. 

The list of solvable matrix models is not very long. Solvability is usally connected to 
reducing the number of relevant degrees of freedom from to N, for matrices of 

size N X N. Two basic and very important examples are given by the one matrix model 
with the partition function: 


ZiMM(a) 


J dM expiV 


1 o (X A 
-tr M2 + - tr 

2 4 


( 1 . 1 ) 


and the two matrix model: 


Z2MM(ct, /3) 


dA dB exp N 


-htvA^ + tv + j{tv A^ + tv B'^) + ^tv{AB) 

^ 4 ^ 


( 1 . 2 ) 

where the integrals go over the N x N hermitean matrices M, A, B. We consider here the 
frequently used quartic non-linearities in the matrix potentials, although one could have a 
more general potential. 

The large N limit of these models is perturbatively equivalent to the planar (or spher¬ 
ical) Feynman graph expansion and describes interesting physical systems. The hrst one is 
related to the problems of enumeration of planar graphs [^] and was proposed as a model 
of pure 2D gravity |^,0; the second, solved in the planar limit in Q and 0 describes the 
Ising model on random dynamical planar graphs p] . It is clear that the exact solvability of 
any natural generalizations or modihcations of these models could be very useful in various 
physical and mathematical applications. 

We present in this paper a very natural solvable generalization of the two matrix 
model given by the following integral: 


Z(a,/3) 


dA dB exp N 



(1.3) 


The only difference with the conventional two-matrix model (|1.2D resides in the last term 
in the matrix potential. This “little” modihcation will turn out to be very important. 

First of all, this new model cannot be solved by the same method as the old matrix 
models: the Itzykson-Zuber-Harish Chandra 0,^ formula is not applicable here to reduce 
it to an eigenvalue problem. In order to avoid this problem, we shall apply a character 


1 









expansion method, which was worked out in connection with the matrix models in the 
papers [§ , ||] , |jT^ where it was successfully applied to the investigation of 2D B? gravity 
@ and the enumeration of branched coverings 0. The basic point of this method is the 
reduction of the degrees of freedom (from ~ to ~ A^) in terms of highest weights of 
representations in the character expansion. It is a sort of Fourier transform for the matrix 
variables. The large N limit allows to apply the saddle point approximation for the highest 
weight distribution (one looks for the most probable Young tableau). Here, following the 
idea of |^, we shall use a double saddle point on both eigenvalues and highest weights. 

Second, the planar graph expansion for the model (|1.3|) corresponds to a different 
statistical mechanical system than the Ising spins on a random dynamical lattice of the 
conventional two matrix model ( |1.2|) . The (/)‘^-type planar diagrams of the model consist 
from the intersecting and self-intersecting closed paths of 2 different colours (see hg. 1). 
The colouring contributes to the entropy and thus describes some spin degrees of freedom. 



Fig. 1: A typical graph in the expansion of the two matrix model with ABAB 
interaction. The two matrices {A, B) are represented by two colours (green, red). 


Alternatively, this statistical model model can be mapped on a special case of the 
8-vertex model on random dynamical graphs. To see it let us set X = A -\- iB-, X is an 
arbitrary complex matrix, and the model can be recast as a one-matrix model: 


Z(6, c, d) 


dXdX"^ exp N 


- J tr(YY^) -f htT{X‘^X^‘^) + ^ tr(A:Yt)2 ^ d ^^^^4 ^ ^^4^ 


(1.4) 

Since X is complex, the propagators carry arrows, and the three vertices reproduce the 
well-known conhgurations of the 8-vertex model. Here, the three corresponding weights 6, 
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c, d are not independent, since they are related to a and jd byS: 

oi fd 


b = 


c — d = 


a — (d 


On a regnlar lattice, the 8-vertex model is critical when the weight d goes to zero, so 
that we are left with the 6-vertex model. In the standard parametrization of the 6-vertex 
model (see for example ||T^ ) 


a = sin — (tt — 6) 

6 = sin ^(tt -I- 9) 
c = sin Att 
d = 0 

where a and b are the weights of two vertices which are indistingnishable after conpling 
to gravity (they are related to each other by a rotation of 7r/4); therefore 9 = 0. For any 
valne of 0 < A < 1, the model has a continunm limit which is a conformal held theory with 
central charg e c = 1, and can be described by a free massless compactihed boson, with 
radins R = \J (np to some dnality transformationsi). 

We expect a similar critical behaviour in our two matrix model, but modihed by the 
presence of 2D gravity [|T^. This is indeed what we shall hnd at the critical point of our 
model, which occurs for = /d^ = l/(47r). It corresponds to the 6-vertex model at the 
point where the weight of (WAfl)^ vanishes: c = 0, or A = 1. The 6-vertex model coupled 
to gravity has already been studied in 0 and shown to be equivalent to the standard 
formulation of the vortex-free compactihed boson coupled to gravity H , with a radius of 
compactihcation 

1 - A/2 


R = 


V2 


which is diherent from the one on the hat square latticei. For our critical point, A = 1 and 
we hnd a radius R = l/(2v^), i.e. one half of the Kosterlitz-Thouless radius. But the 


^ To have the three independent constants, one wonld have to introduce an additional inter¬ 
action 

^ Here, we assume that going from the 6-vertex model to the 8-vertex model correponds to 
perturbing the free boson with a magnetic operator of charge ±4. The Kosterlitz-Thouless phase 
transition occurs at A = 0, as it should be. 

^ This is a sign that the radius of compactihcation is a function of both the weights (i.e. of the 
parameter A) and the geometry of the lattice. We thank I. Rostov for pointing out the discrepancy 
of radii to us. 
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ABAB model is not limited to the 6-vertex model: it also displays off-critical behaviour 
correponding to a particular “slice” of the full 8-vertex model. 

The plan of the article is as follows: in section 2 we shall apply the character expansion 
method to reduce the model to the sum over highest weights of the GL{N) character 
expansion. We then derive the saddle point equations for the density of highest weights 
and eigenvalues, which yield equations dehning the characters and the loop averages in the 
planar limit. 

In section 3 we solve the saddle point equations in terms of incomplete elliptic inte¬ 
grals and hx the parameters of the solution by matching it to the large highest weight 
asymptotics. 

In section 4 we hnd the critical line (of the thermodynamical limit of large lattices) and 
the critical point of the phase transition. We also describe the critical behaviour around 
the critical point, analyse some correlation functions and hnd the string susceptibility. 

The section 5 is devoted to conclusions. 

Finally, the appendices contain a detailed study of three particular lines in our two- 
parameter {a, (3) family of models; they are the lines given by a = 0, /? = 0, and a ~ {3. 


2 . The model and its representation by a character expansion. 

3. L Definition of the model. 


The two-matrix model with ABAB term is given by the partition function 


Z(o,/3) 


dA dB exp N 



( 2 . 1 ) 


where A and B are hermitean N x N matrices. Here a and f3 are positive constants to 
have positive weights in the diagrammatic expansion; of course, as usual in matrix models 
in connection with quantum gravity, the analytic continuation from a, /? < 0 to o, /? > 0 
which allows to dehne the integral o only makes sense at N ^ oo. 

As already mentioned in the introduction, there is no obvious way to do the integration 
over the relative “angle” O between the two matrices A and B, because no formula is 
known for the integral over the unitary group f dO exp[ctr(AOi?Of)^]. To circumvent this 
problem, we use a character expansion of the term exp[A^^ tr(Ai?)^] as a class-function 
of AB\ representations of GL{N) are parametrized by their shifted highest weights hi = 
mi + N — i {i = 1... N), where the are the standard highest weights. Then one has 


{h} 


- l)/2) 

n.('*r"/2)!ni((ft?‘“‘-i)/2)! 


dAdB exp A 


■^(trA^ -FtrH^) -F ^(tr -F trH^) 


X{h}{AB) 


4 







where the sum is over all integers hi that satisfy hi > h 2 >...> Hn > 0, and ^h = 
~ number of boxes of the Young tableau. A(-) is the Van der 

Monde determinant, X{h} is the GL{N) character associated to the set of shifted highest 
weights {/i}, and the ^re the even/odd hi, which must be in equal numbers. It is 

now possible, using character orthogonality relations, to integrate over the relative angle 
between A and B] this leads to a separation into one-matrix integrals: 


{h} 

where is a coefficient: 

1 

and Ri^h^{a) is the one-matrix integral 


( 2 . 2 ) 


R{h}{a) = / dM X{h}{M) exp N 


■-trM^ + -trM^ 
2 4 


(2.3) 


which appears squared in (|2.2| ) because the contributions from the two matrices A and B 
are identical. 


2.2. Study of (a) 

The one-matrix integral R{h} closely resembles what was studied in (see in par¬ 
ticular appendix 2), and we use the same approach. First we go over to eigenvalues: 


R{h}{oi) 


dXk A(A) det j exp N 

k 




(2.4) 


where A(A) = det (^A^ 

Since we now have N degrees of freedom and an action of order we can use a 
saddle point method on the eigenvalues A^. In order to do so, we dehne the resolvent 


w(A) 


tr 1 
iVA-M 



1 

A — Afc 


In the N —>■ oo limit, the A^ form a continuous density on some support [—Ac, Ac], which 
implies that w(A) has a cut on [—Ac, AJ. If we introduce the notation (/’(A) = ^(w(A -1- 
fO) -F(x;(A-fO)), then(/)(Afc) = ^^logA(A). 
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Next we want to differentiate det(A|^^); as det(A|^^)/A(A) is a regnlar fnnction of the 
Afc (since for example, it is a Itzykson-Znber type integral), we can introdnce another 
fnnction h(A) which has the same cut as U!{X) on [—Ac, +Ac] and snch that 

1 o 

H^k) = ( 2 - 5 ) 

(we have introdnced the factor A^ in front of the for convenience). The saddle point 
eqnation of (|2.3|) can now be written as 

A(,^(A)+^(A) - A^+ aA^ = 0 Ae[-Ac,Ac] (2.6) 


Here comes the key remark: as c<;(A) and h(A) have the same cut on [—Ac, Ac], the slashes 


can be removed in (2.6): 

Aw(A) + h^(A)-A2 + aA'^ = 0 (2.7) 

where is h on the other side of the cut [—Ac, AJ. Expanding in powers of 1/A as A —> cx), 
we hnd: 

AT 1 / tr _ \ 

( 2 . 8 ) 


ht(A) = -aA-‘ + V-l-5;^(|A/^" 


n=l 


2.3. Functional inversion and analytic structure of \{h). 

In the same way as we have considered a saddle point on the eigenvalues, we shall use 
a saddle point on the highest weights. Before doing so, we need to understand better the 
analytic structure of the functions involved. 

In a very similar fashion as we have dehned c<;(A), we dehne H{h) to be 

k 


(here, since the hk scale as N in the large N limit, we do not need a 1/N factor in front). 
After appropriate rescaling h h/N, the hk tend to a continuous density p{h) as A —> cx). 
We shall hnd out that part of the density is saturated at its maximum value 1 (the same 
phenomenon occurs in e.g. |]^); therefore we dehne the end points hi and ^2 such that 


p{h) = 1 0 < h < hi 

0 < p{h) <1 hi < h <h2 


Then i7(h ± iO) = i${h) F i'Kp{h) on the cut [0, ^ 2 ], with fi{hk) = log A(h). 


6 






Next we introduce the function L{h) which has the same cut as H{h), and such that 




(the factor of 2 is due to parity reasons). If we dehne X{h) = exp(iL(h)), then, as proven 


in a similar context in the appendix 1 of [0, h(A) and A(h) are functional inverses of each 
other as multi-valued functions. In particuler, on all sheets of A(h) such that X{h) —> cx) 
as h —i> oo, one has, according to 


h = —aX‘^{h) + X^{h) — 1 — 


1 


A2^(h) \N 


tr 


M 


2n 


(2.9) 


A similar expansion was found for large N characters in (see also 0). 

Inverting the expansion (|2.9|) to express A^ (we shall from now on always use A^ and 
not A for parity reasons) as a function of h, shows that there are two sheets that satisfy 
A^ —> cxD as h —> oo: we shall call these A^. One of them, A+ is the “physical sheet” i.e. 
the original function X{h). The simplest analytic structure for A^(/i) is then the following: 
there is a semi-inhnite cut [ha,-t-cx)] connecting A^ and A?_, the hnite cut [hi,h 2 ], and a 
possible pole/zero at h = 0. One can then show that A^ satishes: 

Xl{h)Xl{h) = (2.10) 
a 

The hgure 2 describes this analytic structure. 



Fig. 2: Analytic structure of X^{h). X\{h) has a double zero at h = 0. 
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2.4- Saddle point equation on the highest weights 


We finally consider the saddle point eqnation on {h} in eqnation (|2.2| ). From (|2.4| ) 
one infers that 2{d/dh) logRh = d4{h), while ci^h} has the following large N limit: 


logc^;,} = - y (log y “ ^ log A(h) 

7 / 


SO that, after rescaling of the hi, (d/dh) log(c{^}/A^'^^/^) = —{Ifii{h) + log |)/2. Therefore 
the saddle point equation reads: 


W - ^ ^ log(A//3) Vh€lA,,A2] 


( 2 . 11 ) 


3. Solution of the saddle point equations in terms of elliptic functions. 

In order to solve the saddle point equations, we need to remove the logarithmic cuts of 
the function 2L{h) —H{h) on the physical sheet. This is achieved by dehning the function 
D{h) = 2L{h) — H{h) — 31ogh + log(/i — hi). D has only two (square root type) cuts: 
[hi, h 2 ] and [hs, +cxo], and it satishes the following equations: 


^ = log ^ VhG[hi,h2] 
Ip = log ^ VhG[h3,+cx)] 


(3.1) 


The solution of ( |3.1|) is given in terms of elliptic functions (for some useful formulae on 
elliptic functions, see |T^). Dehne r(h) = \/{h — hi)(h — h 2 )(h — ha) and the incomplete 
elliptic integral of the third kind (h real): 





dh' 

(h-h')r(h') 


(3.2) 


Then D{h) is given by 

D(h) = log^^ - 4^3 (A) + 21>o(A) (3.3) 

Next dehne elliptic parametrizations of h: x is such that 


sn(x, k) = snx — 
8 




hs — hi 
hs - h 













where k = y hl-hl • ^ chosen such that 0 < Rex < X, —iK' < Imx < iK' (where K 
and K' are the quarter-periods), and x{h) has a cut on [hi, -|-cx)]. Also dehne the rescaled 
0 functions: Oa{x) = 9a{TTXj2K), where the 9a, a = 1...4, are the usual 9 functions. 
Finally, let Za{x) = (d/dx) log0a(x). 

Then $^(h) is given explicitly by: 


$^(h) = X 


2Zi(x) 


irr 

K 


+ log 


0 l(x — x) 
0 l(x -t- x) 


(3.4) 


where x = x{h = h). 

Next we must adjust the behaviour of D{h) as h ^ cx), which will hx the unknown 
constants hi, h 2 , hs. Starting from h = —aX'^ -|- — 1 -|- 0(1/A^) (Eq. (PD) , one proves 

that 

2 L±(h) = log — T + ]- + 0{h-^/^) (3.5) 

ct y—ah h 

with L = L_|_. As h —> CX), X ~ \/h^ — hi/A/—h; when x —0, Z’i(x) ~ 1/x, so there is 
a divergent part in (|3.3|) which must cancel out; this leads to the hrst condition: (x(h = 
0) = xo, x(h = hi) = A, x(h = h 2 + fO) = A -|- f A', x(h = h^-\- zO) = iK') 


condition ^ 1: 2xo = K + (3.6) 

TT 

Using ( ^^) and the explicit expressions of and $^ 3 , D{h) can be rewritten in the 
simpler form 


D{h) = log 


h — hi 
—ah'^ 




or alternatively as a function of y = x 

h — hi 


A 

A: 


0 i(xo + x(h)) 


D{h) = log 


-Qch^ 


+ 2 log 


A 


^ 2 {xo+ y{h)) 


(3.7) 


(3.8) 


By matching the Ij^/h terms in the h ^ 00 asymptotics, one obtains the second 
condition: 

condition ^2: Oi = (3.9) 

^ya{hs - hi) 

We have dehned: 

a=1...4 

K 

From (|3.8|) we obtain the expression for the density of highest weights: 


p(h) 


1 + 


log(/3/cr) 

inK 


y{h) - — log 
m 


O 2 (xo - y(h)) 
O 2 (xo+ y(h)) 


(3.10) 


When he [hi, h 2 ], y is purely imaginary, Imy E [0, A'], so that p{h) is real. 
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Fig. 3: Density p{h), q = 0.1, a = (3 = 0.062674. 


To fix hi, /i 2 and h^, we need a third condition, which is the normalisation of the 
density p{h): 


rh^ 

'hi 


dh p{h) = 1 — hi 
Integrating by parts and changing variables from h to y yields: 

Tog(/5/Q;) 


iK' 


dy h{y) 


K 


2{Z2{y - xo) - Z2{y + xq)) 


= —ITT 


(3.11) 


As Z 2 {y — xo) — Z 2 {y + xq) is an elliptic fnnction: 
1 


[Z2iy-xo) - Z2{y + xo)] = -Z4{xq) +cnxo dnxosnxo 


dc^y 


1 — sn^ Xo dc^ y 


(dc = ^). one can express ( |3.11 ) in terms of the complete elliptic integrals of the hrst 
and second kinds. The resnlting eqnation hxes a: 


condition ^ 3: 


a 


^ 'e'Qi - K'Qi ^ 




sn^ Xo 


(3.12) 


To snmmarize, given \og{(3/a) and the elliptic nome q = e , one can compute xo 

from (|3.6|) , then a and /3 from ( p.l2|) , and hnally hi, h 2 , h^ from ([ 


4. The critical line and the critical point. 

We now investigate the critical properties of the model. Just as in the standard 
eigenvalue problem, criticality is attained when a branch point collides with a critical 
point (i.e. a point where the derivative vanishes), so that the square root singularity 
degenerates into a (h — ho)^/^ behavior. If one starts with small coupling constants a and 
fd and then increases these constants, two phenomena can occur, which correpond to two 
pieces of the critical line. 
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4 . 1 . Criticality of type A. 

The first type of criticality appears when the singularity appears in the vicinity of ^ 3 , 
i.e. at the start of the inhnite cut [h 3 ,+cx)]. This is what happens, for example, if one 
hxes /? to a small value, and increases a: for a smaller than a critical value ac(/?), is 
well-dehned and real, but for a > Uc becomes complex, which is the sign of a change 
of analytic structure, and therefore of criticality. 

In order to hnd the critical line in the {a, (3) plane, one looks at the behavior of D{h) as 
h ^ h^, i.e. X —> IsiiK'. The procedure is very similar to looking at the a; —> 0 behaviour; 
criticality is obtained by cancelling the \/h — term, that is when 


Sl4=!^* + 4Z4(xo)=0 
A 


(compare with (|3.9|) ). 

Using (|3.9[) one can simplify this condition: 


cna:o dna:o 


snio i\/a{h3 - hi) 


(4.1) 


and reexpress it in terms of h: 

hih2 1 

^3 16a 


(4.2) 


4 . 2 . Criticality of type B. 

This time the singularity occurs at the level of the hnite cut [hi,h 2 \\ this is what 
happens for small a and increasing /?. 

Criticality is more precisely found by looking at the behaviour of D{h) as h —>■ h 2 , i.e. 
X ^ K ± iK'. Cancellation of the y/h — ^2 term implies: 

ll3 = i^* + 4Z3(x„)=0 (4.3) 

Combining it with ( p.9|) yields 

cnxo 1 

sn xo dn xq 4 ^Q;(h 3 - hi) 


or equivalently 


hih^ _ 1 

^2 IbCK 


(4.4) 
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4-3. Double criticality: the critical point. 


When one imposes both ([4.21) and (|4.4|) , one hnds = ^ 2 , hi = l/(16a): the critical 
point correponds to the limit when the two cnts tonch each other. Therefore we are in the 
trigonometric limit /c —> 0. We also hnd a = j3. In this case one can solve directly the 
eqnations (|3.1|) which combine into a single eqnation; the resnlt is: 


D{h) = log -—^ - 2 log = log-—— -41og(Vhi - h + ^/hl) 

—ah^ —y/hi — h + y/hi —ot 


p{h) 


1-arctan \/h/hi — 1 

TT 


p{h) is zero for h 
density gives 


2hi, so one has h 2 = hs = 2hi. Finally the normalisation of the 



and therefore a — jS — l/(47r). 

The critical line and the critical point have been drawn nnmerically on hgnre 4. The 
critical line is divided in two: the left of the critical point (a > /3) is the phase A (criticality 
is of type A), whereas the right (/3 > a) is the phase B. Note that the trigonometric limit 
g —0 is precisely the neighborhood of the critical point, and on the contrary that the 
hyperbolic limit q ^ 1 Is the limit a —0 or 0. 



0.00 0.05 0.10 0.15 0.20 0.25 

P 

Fig. 4: Critical line, critical point and equipotentials of q {q = 0(0.03)0.6) in the 
{a, (3) plane. 
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Fig. 5 : hi, h 2 , hs on the critical line as a function of (5. 


To complete the picture at the critical point, we need to go back from D{h) to L{h), 
which requires to compute H{h). In fact, H{h) is not a standard function but dH/dh is: 


dH 

Ih 



p{h') 


dh' 

[h-h'Y 


2 1 I hi ^ y/h — hi + y/hi 
Trhyh — hi ^ y/h — hi — y/hl 


This expression is sufficient to expand 2L{h) = D{h) — i7(h) + 3log h — log(h — hi) around 
the “critical point” he = h 2 = h^, which provides us with the loop function h(A) for large 
loops. Only H{h) has a singularity at this point: dH/dh oc log(h — he), and therefore 
H{h) (X {h — he) log(h — he), or hnally 


A — Ac oc (h — he) log(h — he) 


which we can invert: 


h — he X 


A-Ac 

log(A - Ac) 


(4.5) 


We shall now study the properties of the vicinity of the critical line. In particular, we will 
need the exact expression for a correlation function. By functional inversion, A(h) provides 
us with all correlation functions of the type (tr In practice, one can easily compute 

(tr by inverting (|2.9|) up to order 1/A^: 


2L(h) = log —- — 
a 


1 1 

y/—ah h 




+ 


24^3/2 2^J 2(-h)3/2 


+ 0(h-") 
(4.6) 
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where L{h) = 21ogA(/i), and expanding the explicit expression (|3.7|) of D{h). 
Here, we shall nse another correlation fnnction, (tr It is given by: 




= —Z 
d(3 


1 

1 


1 


-hl-- + 
2 ^ 2 


1 

2 


nK' 


dy 


p{h)hdh 


dp 


dy 2 


(4.7) 


At this point, the calculation becomes similar to that of a. The result is: 


?»»'■> " 7, 


+ hs + TT- 


E'(04/sn^ xo - 04(2 - A:'")) + A:' 04(1 - k'^) 


(E'Oi - A:'04/sn2 a:o)2 


(4.8) 


4.4- Vicinity of the critical line. 

Let us suppose hrst that we approach the critical line at hxed slope s ^ Q\ then q 
tends to a non-zero critical value qc, and one can show the following facts: -^oi\q=q^ = 0, 
^ <^tr(Ai?)^)^_^ = 0 and also -^hk\q=q^ = 0 with /c = 1, 2 for criticality of type A, /c = 1, 3 
for criticality of type B. This implies in particular that 

-{abA 

^ / sing 

where A = ckc — a is the renormalized cosmological constant, and therefore the string 
susceptibility exponent is — this means that on the whole critical line except at the 
critical point, the continuum theory is a c = 0 theory (pure gravity). 

In fact, the analytic structure found around any point on the critical line in phase A 
(resp. B) is identical to that found for the endpoint a = 1/12, /3 = 0 (resp. ft — 2/9, a = 0). 
These two particular points have been studied in detail in appendix A. The conclusion one 
draws from the more explicit calculations done at these points is the following: though 
we have in both phases a c = 0 theory, for phase A the loop scaling f un ction is the usual 
pure gravity loop function, whereas in phase B this function describes non-trivial loops 
(eq. (A.6)). 
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4-5. Vicinity of the critical point. 

We now investigate the vicinity of the critical point a — (3 — l/(47r), which is also the 
small q region. We refer to appendix B for details on small q expansions. 

If we dehne the deviation from criticality 

A = 1 — dTTQ! 


and the slope in the (cr,/?) plane 

^ log(/?/cr) 

47r 

then we obtain up to order q^ the equation dehning q in terms of a and (3: 

A = 2s + log q~"^ + 2s^ log^ + 4^^ — 4s^ (4.9) 


We can hnd from here the shape of the critical line near the critical point. It is given 
by the equation: 

d 4s‘^ 

—A = 16q\ogq ^ -logg ^ = 0 (4-10) 

dq q 

or 

q = ±s/2 (4.11) 

for /3 < a and (3 > a, respectively. One can check (using the last two expansions of (B.3)) 
that it hts well the conditions ( [I.l|) and (|4.3|) . 

In particular, this implies that the slope of this line at the critical point is: 


da 1 
dp 1 — 271 


(4.12) 


It is possible to expand (^T^) around Q = 0; the result is of the form 


tr 

— {AB)‘^ ) =Ci + C2S + C3S^log^g“^ +C4S^logg“^ +C5q‘^\ogq~^ + CqS^ + Crq'^ (4.13) 


where the coefficients Cj are given in appendix B. 

For s = 0, from ([4.9|) A cx logg“^, and the singular part of ([I.13|) yields 


oc ^ 


(4.14) 


that is a zero string susceptibility exponent, plus logarithmic corrections. This is charac¬ 
teristic of a c = 1 model coupled to gravity |T^ . 
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Let us also note that on the critical line s = ±2(j, one can check that, as expected, 


dq 



= 0 


The free energy of the underlying statistical-mechanical model in the thermodynamical 
limit can be obtained by the following trick 0 . The planar free energy of the matrix model 
(i.e. the partition function of the statistical model on connected planar graphs) near the 
critical line looks as 

logZ ~ (4-15) 


where q:*(s) is the radius of convergency of the series at hxed s, and n (the number of A'^ 
vertices) characterizes the size of the graphs. Therefore, log[a*(s)] plays the role of the 
free energy per volume unit in the thermodynamical limit. By plugging ( [4.11|) into (|4.9| ) 
we hnd: 

A*(s) = 1 — 47rQ!*(s) = 2s log(4/s^) s^ log^(4/s^) — 3s^ (4.16) 


Whereas this function and its hrst derivative are hnite and continuous the second 
derivative (playing the role of the specihc heat of the problem with s along the critical line 
being the “temperature”) developes a log^ type singularity: 


ds'^ 


A*(s) 


~ 21og^(4/s^) 


for 


(4.17) 


giving rise to a second order phase transition. 

This is an unusual result since for most of the spin models on random dynamical 
graphs the transition is usually of the third order. A similar situation, but with a log type 
singularity accurs in the model of percolation on dynamical random graphs solved in . 


5. Conclusion 

Our results show that the possibilities to hnd new exact solutions of the matrix models 
are far from being exhausted. The character expansion method used here enlarges the class 
of such models with a potentially very rich spectrum of critical behaviours and thus with 
numerous possible applications to mesoscopic physics, enumeration of graphs, statistical- 
mechanical models on random dynamical graphs, etc. 

For example, a rather general class of such models for p hermitean matrices can be 
given by the action: S = tr(^^^j^ Vk{Mk) + f{Y[k=i ^k)) for arbitrary functions 14 and 
F. It is easy to see from the orthogonality of characters that the exponent of the last 
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term can be integrated over relative angles of the matrices. But even this is not the most 
general example. 

One interesting question we can pose here is the following: is there any integrable 
structure behind the representation expansions of partition functions like (|2.2|) ? For more 
conventional matrix models, we are used to the fact that they can be represented as 
tau-functiopns of some integrable hierarchies of classical differential equations. These tan- 
functions are usually some N x N determinants or Fredholm determinants in the grand 
canonical ensemble with respect to N. 

The particular elements of this formula remind the objects known from the tan- 


functions: the characters or their integrals with one-matrix weights (like in (|2.3D ), as well 
as the product have determinant representations (or, in other words, 

some fermionic analogues). But we don’t know such a representation for the whole sum 
over Young tableaux. If such a representation exists it could allow us to analyse the critical 
behaviour for all genera of random graphs in the double scaling limit, like for the more 

1 , 0 . 


standard matrix models 21 


Let us also note that the universal behaviours of the physical quantities we found, 
like the loop function (|4.5|) look different from the known patterns for the c = 1 matter 
coupled to 2D gravity. A possible explanation is that the loops we are considering differ 
from the usual ones (through different “boundary conditions” on the loop). If so, it would 
be interesting to make the connection between the different types of loops. 

Also, the critical point we have in our model is only one point in the critical line of 
c = 1 CFTs. In order to describe more precisely the properties of this point, it would be 
interesting to investigate the possibility of a non-trivial scaling limit around the critical 
point CK* = /3* = 1/ (47r). This would shed light on the perturbing operator around the hxed 
point, which is directly related to the radius of compactihcation R in the equivalent bosonic 
formulation; it would give an independent check of the predicted radius R = 1/(2a/2). 

Finally, it would of course be very interesting to obtain the solution of the full 8-vertex 
model coupled to 2D gravity (with all couplings b, c, d independent in ( |L^) ). This would 
allow to describe the most general critical properties of c = 1 coupled to 2D gravity. But 
for the moment, even character expansion methods cannot solve this more general matrix 
model. 
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Appendix A. Study of particular cases. 

A.l. /3 = 0; the usual one matrix model revisited. 

When /3 = 0 the model reduces to two decoupled one-matrix models with the standard 
action i trM^ -|- ^ tr We shall now compare our solution with the classical solution 
of pure gravity. 

Since the expansion in characters is trivial for (3 = 0, the saddle point character is 
necessarily the trivial representation, i.e. p{h) = 1 for h G [0,1]. Therefore 

H{h) = 


which implies that hi = ^2 = 1 (the cut reduces to a pole at h = 1), and 


xl{h) = 1 

a/o y/hs - 1 =F y/hs - h 
hs is given explicitly as a function of a by 


(A.l) 


, 5/4 - 3a + y/1 - 12a 

AlQ = - 

9a 

which has a singularity at a = 1/12, the critical point of pure gravity. 

On the other hand, one can use the standard method for this model, which is to 
solve directly the saddle point equations for the eigenvalues. Indeed, when the character 
is trivial, it is clear from (|2.5| ) that h(A) = Aci;(A); then ( |2.6|) becomes 


2^ - A^ + aA^ = 0 


which yields 


h± 


A2 - aA^ ± (aA^ - 1 + l^aXl)y/ 


(A.2) 


where h(A) = h_|_(A) is the physical sheet, and {X) = h_(A) is the other sheet connected 
by the cut [—Ac, +Ac]0. The end point of the cut Ac is given by: 


A 


2 

c 


2 

3a 


(1 - y/1 - 12a) 


Note that {X) has the correct expansion (|2.8|) as A —> cxd. 


^ As a function of A^, the cut of h{X^) is [0, A^]. 
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It is easy to show that A^(h) given by ( |A.1|) , and h(A^) given by (|A.2|) , are fnnctional 
inverses of each other: indeed they are the solntion of the eqnation 


aX^{h - 1 ) - X^{h - ahsXl) + h‘^ = 0 


(A.3) 


which is quadratic in A^ and in h. 

Criticality is of course obtained when a —> cTc = 1/12. By computing correlation 
functions one immediately hnds that the string susceptiblity 7 = — that is pure gravity. 

Let us also remind the reader of the universal loop scaling function. The asymptotics 
of the loop average (Iy (which is interpreted as the summation over surfaces with 
a hxed boundary of length n), n large, are dominated by the singularity of the resolvent 
c<;(A) (or of h(A)) which is closest to A = cx). Here, it is the square root singularity starting 
at A = Ac. Equivalently, one can say that he = h{X = Ac) is a critical point of the function 
X{h). When a cTc, he collides with the branch point ha, as expected. One can now take 
a scaling limit where the renormalized cosmological constant A = ae — a 0 (so that the 
average area of the closed surface (A.) oc 1/A diverges) and the renormalized boundary 
cosmological constant A — Ac —»• 0 (so that the typical size of the loop (/) oc 1/(A — Ac) 
diverges) in a correlated ma nn er. Here, the scaling function is: 


h he f X Ac \ f ^ Ac 
A3/4 “ Al/2 ) AV2 


1/2 


(A.4) 


A.2. a = 0. 


When q; = 0, only the vertex tr{ABAB) survives, so that one has a model with two 
types of loops which intersect each other all over the surface. Discarding one of the two 
types of loops, one hnds a model of dense loops, i.e. loops that cover the entire space. 
This looks very similar to the dense phase of the 0(1) model (i.e. low-temperature phase 
of the Ising model coupled to gravity), which has already been studied in [^. In fact. 


the ABAB model with a = 0 (and its 0(n) generalisation) was already studied in 


(by completely different methods from ours), where it was shown to be equivalent to a 
modihed 0(1) model. Let us compare here our results with those of and . 

If q; = 0, the expansion ( |2.9| ) implies that the analytic structure looks different from the 
generic case cr > 0: indeed, A has no cut going to inhnity. In other words, when a ^ 0, ha 
tends to inhnity (more precisely, one can show that h^ ~ 1/(4q;)). We are left with a single 
cut [/ii,h 2 ]; using the general solution and the additional relation A^(h) = hexp{H{h)) 
which follows from the fact that a = 0, we can explicitly the function A(h): 


X{h) = 


1 {h y/{h — hi)(h — /i2) + \/h]h^)‘ 

2 h{\/h — hi + y/h — ^ 2 ) 


(A.5) 
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hi and /i 2 are given as solntions of a third degree equation: hi ^2 = “*■ ~ I where 


X 


16 64 

:X + = 0 


3/32 3/32 


The discriminant A' = — 1 vanishes when (3 — f3c — ^ (the value of f3c coincides with 

what was found in [^). Then hi = |, ^2 = and Ac = \{h = /i 2 ) = 3/a/2 satishes 
PAl = 1 , as follows from general arguments ||25|| . 

One can easily compute the string susceptibility in the limit P ^ Pc', using x = 
6(1 + -^A '^/2 _|-one hnds 


/ SI 


(X 


A'3/2 


Sing 


so that the string susceptibility exponent 7 = — 1 / 2 , which indicates a c = 0 model. 

However, the loop scaling function is non-trivial; in order to see this, we use the simple 
formula: 

^ h - 2y/hih2 

—Lih) = 2A — =- , 

dh ^ ' dh h^/{h- hi){h- h2) 

which shows that X{h) has a critical point at he = 2^/hlh2. This point meets the branch 
point h 2 when P —> Pc, as expected. In the limit P —> Pc, the proper scaling ansatz is that 
h — he (X A'^/ 2 ^ A — Ac oc A'^A. Then the universal loop scaling function is given by: 


A-Ac 

A'3/4 


( h — he 
[ a'V 2 


2 


^cst + 


h — hc\ 
A'1/2 J 


1/2 


(A.6) 


Note that the scaling A — Ac oc is unusual since it implies that (/) oc 

In both limits ^ 7 i /2 —^ 0 and ^ 7 i /2 —*■ 00 ? oue recovers the same asymptotics as the 
loop function found in 1 ^^; however, the full loop function is different, suggesting that for 
loops of the same size as a typical loop on the surface, the physics of these two models 
differs qualitatively. 


A.3. a = p. 

This is the line on which the critical point lies; on this line it is possible to reexpress 
the function D{h) using only standard functions (log, square root ...). However, the 
parameters of the problem still depend on the constants k, K', E' which are all functions 
of the elliptic nome q. Let us see how this works. 
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By putting a ~ (3 we obtain from (|3.6|) that xq = Kl2 and from (^) that 

k' 


hi = 

h2 = 

h, = 


4q;(1 + k') 
k' 

4a(l + k') 

1 

4q:(1 + k'Y 


(A.7) 


If we employ the functional relations for sn and cn, together with the formulas siv{K/2) = 
7 ^ etc, we hnd: 


2 {h - hi){h2 -h) + ^(^2 - hi){hs - h) 

pin) = 1 + — log — , - , 

*7r [h — hi){h^ — h) + {h^ — hi) (^2 — h) 

The constants hi, h 2 , hs are given by (|A.7|) , where a is: 

E' -{1- k')K' 


(A.8) 


a = 


27r(l + k') 


(A.9) 


Appendix B. Small q expansions. 

We rewrite the condition (|3.6|) as 


TT TT , _i 

—xo = j+slogq 


where it is reminded that 


s = 


log{/3/a) 


47r 


In these variables, the dehning equation ( |3.12|) for a becomes: 


27r^cr = KE'- 


1 


e'. 


2c -u 'i 

^ oA 2 K ) 


-KK'- 


2 « + it(if) 


2^+77 (it) 


We then use the q expansions: 


K = Ti/2{l + Aq + Aq^+ 0{q^)) 

K' = —(1 + Aq + Aq^ + 0(i?^)) 

E' = 1 A- Aq log q~^ — Aq — 8g^ log q~^ + 12q^ + 0{q^ log 
k = Aq^^‘^{l — Aq + lAq^ + 0{q^)) 
k' = l-8q + 32^2 ^ 


(B.l) 


(B.2) 
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and the double q, s expansions: 


l/sn^(a:o) = 2(1 — 2slog(j ^ — 4^ + 4s^ log^ g ^ + 16q'slog(j ^ + 16q'^) + ■ ■ ■ 
^ (i?) = 1 “ 2s log + 2s^ log^ q~^ + 4q2 H- 




^ /TTXqX 

03 \2 k ) 


-4q H- 


(B.3) 


where s is assumed to be - at most - of order q: this turns out to be a correct assumption, 
since one has |s| < 2^ (the limiting values dehne the critical line). In ( |B.3|) the error is of 
order q^ log^ q~^. 

The deviation of a from criticality is then given by: 


A = 2s + logg ^+2s^log^(j ^ + 4(j^ — 4s^ + 0((j^ log^ g (B.4) 


The analogous expression for A' = 1 — 47r/3 is: 

A' = 2(1 —27r)s + 8q'^log(j“^ + 2s^log^(j“^+4q'^ —4s^(l —27r + 27r^)+0((j^log^(j“^) (B.5) 

On the critical line, s ~ ±2q and therefore its slope is A/A' 1/(1-2x). 

In the same way, one can expand the correlation f un ction (|4.7|) ; the result is 

^ — 3 + 27r + 2s(—3 + Qtt — 47r^) — 6s^ log^ q~^ + 247rs^ logg”^ 

— 24^^ logg“^ + 47rs^(3 — 127r + 47r^) + 12q'^(—1 + 47r) j + 0{q^ log^ q~^) 

(B.6) 


which is of the form ( 4.13 ). 
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